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states at the (multi-)TeV scale, beyond the direct reach of the LHC. 

Keywords: Higgs Physics, Effective Field Theories 


ArXiv ePrint: 1506.05465 



Contents 


1 Introduction 1 

2 Experimental sensitivities 3 

3 Scalar UV completion 5 

4 Singlet-doublet UV completion 7 

4.1 Effective field theory approach 8 

4.1.1 Precision electroweak observables 13 

4.1.2 Higgsstrahlung 17 

4.2 Loop computation of EWPO 22 

5 Conclusions 23 

A One-loop computation of S and T 25 

B A clarification on matching in the presence of mixing 27 

B.l Example 1: two-point functions of the UV-theory m field and the EFT field y 28 
B.2 Example 2: Z-boson coupling to the EFT field y 29 


1 Introduction 

Among the particularly compelling scenarios in which new physics couples to the Standard 
Model (SM) are the so-called vector, neutrino and Higgs “portals,” which involve the three 
lowest-dimension SM-singlet operators which can couple to new physics (see, e.g., ref. [1] 
and references therein). Of these, the Higgs portal, which encompasses operators of the 
form where Onp is a SM-singlet operator built from new fields, has been the 

focus of much interest in the literature (e.g., refs. [2-15]), especially in the wake of the 
2012 discovery [16, 17] of a light Higgs boson with SM-like properties [18, 19]. There are a 
number of reasons for its appeal aside from its low dimensionality. The Higgs portal provides 
a possible coupling of the dark matter (sector) to the SM, as we discuss in more detail below. 
Some well-studied models of this class—such as H^HS^ with S a scalar singlet—have the 
ability to effect non-trivial modification of the dynamics of electroweak symmetry breaking 
(EWSB), rendering a first-order phase transition possible, with interesting implications 
for electroweak baryogenesis [20-25]. Additionally, such models may characterize the low- 
energy physics in theories which ameliorate the SM hierarchy problem [26]. 

While the simple picture of only considering the lowest-dimension fermionic Higgs por¬ 
tal operator in a bottom-up analysis is of course instructive, much more can be said about 
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Figure 1. Two possible UV completions of the fermionic Higgs portal operator, H^HxX- 
the “scalar completion” (left diagram), a singlet scalar S has a renormalisable coupling to the H 
doublet and a Yukawa coupling to the singlet x held; in the “fermionic completion” (right diagram), 
a Yukawa-like coupling exists between a vector-like SU{2) doublet F and the H and x fields. 


particular possible UV completions of this model. In this paper we study two such com¬ 
pletions (see figure 1) for the CP-even fermionic Higgs portal operator, Hxx^ their 
precision constraints. In the first model, a heavy singlet S acts as a mediator, with a Yukawa 
coupling to the fermion x (which we take to be Dirac) and renormalisable couplings to the 
Higgs-sector of the SM. This is essentially an extension of the singlet-scalar Higgs portal 
scenario. In the second model, we introduce a new heavy vector-like SU (2) fermion doublet 
F and couple both the singlet and doublet fermions to the SM Higgs: —kHFxF^-c.. Such 
singlet-doublet models have also been studied in the past [27-29]. 

Since the new states in either completion are coupled to the Higgs-sector of the SM (and 
in the latter scenario, also have non-trivial charges under the SM gauge group), precision 
electroweak and Higgs physics will be important probes of these potential UV completions; 
for example, measurements of S, T and the Higgsstrahlung cross-section azh are particularly 
sensitive. Although the present electroweak precision constraints have been known for a 
long time to be stringent, the current round of conceptual studies for future precision high- 
energy e^e“ machines—e.g., the ILC, CEPC, and FCC-ee^—have resulted in significantly 
sharper projections for how the constraints could improve. This brings into even clearer 
focus the discovery and exclusion potential possible through a study of these probes. It 
is the major goal of this paper to more fully understand how the projected sensitivities of 
these machines will allow stringent constraints to be placed on the UV completions of the 
Higgs portal that we consider. 

Although we do not impose a dark-matter interpretation on the singlet x in our models, 
it is worth noting this possibility has received much attention recently. For the dimension- 
five fermionic Higgs portal, H^H{axx + which comprises both CP-even and CP- 

odd coupling terms, a thermal-relic dark-matter interpretation is tightly constrained in the 
CP-even scenario by spin-independent direct detection constraints, whereas for the CP-odd 
coupling the direct detection constraints are alleviated and much parameter space remains 

^The collider formerly known as TLEP. 
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open (see, e.g., refs. [8, 15, 30-32] and references therein). Turning to the UV completion of 
these operators, additional parameter space for thermal-relic dark matter is available when 
the new mediator mass (either the doublet fermion or singlet scalar) is similar to or below 
the mass scale of x- For example, the relic density may be set by coannihilation of the 
charged and neutral fermions [29, 33], or by dark matter annihilation into lighter scalars 
(see refs. [7, 34-37] for studies of the singlet scalar UV completion). These scenarios can 
have small couplings and thus evade electroweak precision and other constraints. As our 
focus is not on the dark-matter interpretation of X) but rather on precision probes of the 
fermionic Higgs portal coupling and its UV completions generally, we restrict ourselves to 
the regime where the new mediator is heavy and the EFT is valid. 

The remainder of this paper is structured as follows: we begin in section 2 summariz¬ 
ing the experimental sensitivities on measurements of S, T and azh which are currently 
projected to be attainable at future colliders. In section 3 we introduce and briefly discuss 
the singlet-scalar UV completion, summarizing existing results in the literature. In section 
4 we introduce the fermionic UV completion and discuss our general expectations. We then 
detail the construction of the EFT we used to analyse this model in section 4.1, and give the 
EFT results for the precision electroweak and Higgsstrahlung constraints in sections 4.1.1 
and 4.1.2, respectively. In section 4.2 we discuss the one-loop computation of the precision 
electroweak limits in the general mass parameter space for this model. We conclude in 
section 5. Appendix A contains technical details on the one-loop computation of S and T, 
while Appendix B contains a clarihcation on a technical point of the EFT matching. 


2 Experimental sensitivities 


Several proposed e'’'e“ collider experiments which can make advances in precision elec¬ 
troweak and Higgs physics measurements are currently under consideration: the Inter¬ 
national Linear Collider (ILC) [38], the Future Circular Collider (FCC-ee) [39] and the 
Circular Electron Positron Collider (CEPC) [40]. We summarize the sensitivities for vari¬ 
ous present and future measurements of the electroweak S and T parameters [41, 42], and 
the Higgsstrahlung cross-section azh^ along with the references on which they are based, in 
table 1. For each collider, we include several scenarios for the run parameters in order to 
compare the improvements in sensitivity that are possible with collider upgrades. 

For the S and T electroweak precision obserable (EWPO) limits, we present a parametri- 
sation of the covariance matrix used to construct the 68% coverage likelihood contours pre¬ 
sented graphically in the references. We made the simplifying assumption of exact Gaussian 
uncertainties centered at {S,T) = (0,0), and have parametrised the covariance matrix as 


S = 


cr: 


PsT 


PsT 


( 2 . 1 ) 


from which the likelihood function is given by —21og[T/To] = where A^ = 

(S', T). The parameters (cTs, ctt, Pst) were obtained by least-squares htting of the constraint 
equation —2 log[£/£o] = 2.30 to a large number of (S, T) co-ordinates read from the relevant 
graphically-presented 68% coverage ellipses. 


Table 1. Experimental sensitivities for current and future measurements of {S,T) and azh 
assumed in this paper. The azh sensitivities are quoted as 68% confidence percentage uncertainties, 
while the {S, T) limits are shown as the parameters ( CTs, Uti Pst ) which define the covariance matrix 
in eq. (2.1), and are extracted from the 68% coverage likelihood contours (—21og[£/£o] = 2.30) in 
the indicated references, f The individual ILC Higgsstrahlung constraints must be combined where 
necessary by summing the relevant — 21og[£/£o] values to obtain the constraints for the three 
scenarios we consider in this paper: 250/fb @ 250 GeV, 250/fb @ 250 GeV + 500/fb @ 500 GeV, 
and 1.15/ab @ 250 GeV + 1.6/ab @ 500 GeV. 


Precision Electroweak {S,T) 


Scenario 


( <Ts , <Trp , Pgrp ) 

Reference 

Cnrrent 


(8.62x10-2 , 7.37x10-2 , 0.906) 

[43], hgnre 1 

CEPC “Baseline” 


(2.39x10-2 , 1.93x10-2 , 0.844) 

[43], hgnre 4 

CEPC “Improved Fz, sin^^” 


(1.14x10-2 , 8.79x10-3 , 0.518) 

[43], hgnre 4 

CEPC “Improved F^, sin^0, 

m/’ 

(1.12x10-2 , 7.26x10-3 , 0.779) 

[43], hgnre 4 

ILC 


(1.71x10-2 , 2.14x10-2 , 0.891) 

[43], hgnre 1 

FCC-ee-Z (aka TLEP-Z) 


(9.32x10-3 , 8.70x10-3 , 0.440) 

[43], hgnre 1 

FCC-ee-t (aka TLEP-t) 


(9.24x10-3 , 6.18x10-3 , 0.794) 

[43], hgnre 1 

Higgsstrahlung crzh 

Scenario 


{^crzh)l(^zh (%) 

Reference 

CEPC: 5/ab @ 240 GeV 


0.5% 

[40] 

ILC: 250/fb @ 250 GeV 


2.6% 

[44, 45] 

ILC: 500/fb @ 500 GeV 


3.0% 

[44, 45] 

ILC: 1.15/ab @ 250 GeV 


1.2% 

[44, 45] 

ILC: 1.6/ab @ 500 GeV 


1.7% 

[44, 45] 

FCC-ee: 10/ab @ 240 GeV (4IP) 

0.4% 

[39] 


Tn this parametrisation, the T-axis intercepts of the 68% coverage likelihood contour for two parameters 
are at 5 = 0, T = ±(Tt[1 — Pst] • [(^X^) 2 *^] where = 2.30, and similarly for the S'-axis 

intercepts. In particular, this means that for S' = 0, the coverage in the T parameter is about a factor of 2 
better when going from the FCC-ee-Z to the FCC-ee-t scenario, despite the modest 30% decrease in a^. 


The precision electroweak sensitivity is based primarily on Z-pole measnrements at 
the varions colliders. For the CEPC, the “baseline” is a Z mass threshold scan. The 
“improved” scenarios consider the gains if the collider is npgraded to allow for improved 
measnrements of T^ (with the resonant depolarization method for energy calibration), 
sin^ 6, or possibly rrit (from an ILC top threshold scan). The FCC-ee-Z sensitivities are for 
Z-pole measnrements with polarised beams (denoted “TLEP-Z” in ref. [43]). The FCC-ee-t 
scenario inclndes Z-pole measnrements with polarised beams along with threshold scans 
for WW and tt (denoted “TLEP-t” in ref. [43]) and has improved sensitivity owing to more 
precise measnrements of the W and top masses. 

As has recently been emphasized in refs. [26, 46] (and investigated farther in, e.g., 
refs. [47-51]), the “Higgsstrahlnng” process e^e“ —?• Zh can be measnred at the percent or 
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sub-percent level at future colliders and is a sensitive probe of new physics.^ Among the 
various ways a new-physics model modifies the Higgsstrahlung cross-section, a modification 
to the wavefunction renormalization of the Higgs (i.e., a modification of the momentum- 
dependent part of the Higgs two-point function) may be induced, leading to a “Higgs 
oblique” correction. 

The limits on zh!^Zh are quoted in table 1 as the 68% percentage confidence bounds 
given in the various references. We consider a CEPC projection with 5/ab of data at \/s = 
240 GeV [40] and an FCC-ee projection with the assumption that there are 4 interaction 
points with 10/ab of combined integrated luminosity collected at = 240 GeV [39]. 
Various configurations of proposed runs for the ILG have been analyzed in the literature 
[44, 45, 58]; we consider only the configurations where (a) 250/fb of data are collected at 

= 250 GeV, (b) 500/fb of data at -^/s = 500 GeV are added to the 250/fb of data 
collected at ^fs = 250 GeV, and (c) the integrated luminosity is increased to a total of 
1.15/ab at \/s = 250 GeV and 1.6/ab at \/s = 500 GeV. 


3 Scalar UV completion 


The simplest possible UV completion of the GP-even fermionic Higgs portal is to take 
the Standard Model (SM) augmented by a vector-like Dirac fermion SM-singlet y and a 
SM-singlet scalar S. These couple via the following Lagrangian 

C = CsM + ix$X - m^XX 

+ - ^msS^ - + amsS\H\^ + - KgSxx, 

(3.1) 

where a, 5, e^, ks, ms and m^. are real parameters. To be explicit, our sign and normalisation 
conventions for the SM Higgs-sector are 

+ (3.2) 

and our unitary-gauge normalization conventions are such that H = (0, {v + h)/\/2), with 
V ^ 246 GeV. 

This model is just the renormalisable scalar Higgs portal model (see, e.g., ref. [50] 
for detailed discussion closely related to the present work), augmented with the singlet-y 
Yukawa coupling. If we consider the limit where bms,€s 'C ams,Ks, and further assume 
that ms S> v,m^, the mediator S can be integrated out at tree-level to give rise to the 
following effective field theory: 


Teft a S 


+ 'tX0X - m^XX 

1 


+ + \'\{xxf + • • 


m^ 2 


2 


• , (3.3) 


^See also refs. [52-57] for previous work on constraining anomalous Higgs couplings with the Hig¬ 
gsstrahlung cross-section. 
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where • • • represents terms proportional to various powers of €s,b and neglected higher¬ 
dimensional operators. Allowing terms proportional to €s,bms leads, at dimension-6, to 
only one additional operator, \H\^, which taken together with the ability to change the sign 
of the operator in the SM-Higgs Lagrangian and yet maintain a stable minimum to 
the Higgs potential, can have interesting implications for the order of the electroweak phase 
transition [50, 59]. 

For the purposes of the discussion here, the \H\^ operator can simply be absorbed 
into an unobservable shift of the SM Higgs quartic-coupling: A —)• A -|- The fermionic 
Higgs portal operator H'^Hxx leads to a variety of effects (see, e.g., ref. [31] for a detailed 
analysis): upon EWSB, it contributes to the mass of the x held, and allows both hxx 
h^xX couplings. The former coupling allows for X'X scattering via Higgs exchange, and 
there is a further scattering contribution from the (xx)^ operator. 

It is possible in this scenario to have the x held play the role of the stable dark matter 
(DM), saturating the relic density. If the mediator S is assumed to be somewhat heavier 
than the X) EFT analysis is applicable [30, 31]: the relic abundance requirement fixes 
ms/{aKs) ~ 0.5 TeV, while satisfying the LUX direct detection bounds [60] requires that 
the X be fairly massive, ^ 3 TeV. These requirements are actually in tension with 

the assumptions of an EFT analysis as they demand a very large coupling ans ^ ^vr to 
avoid > ms and the resulting issues with perturbative unitarity [61-63] for the EFT 
description of the non-relativistic DM freeze-out. This essentially means that, with the 
finely-tuned exception of ~ [30, 31], the heavy-mediator scenario is ruled out 

here. On the other hand, the scenario with a light mediator, mg < is still viable as the 
direct detection constraints can be avoided while maintaining the correct relic abundance 
[30]. Of course, if we were to abandon the DM interpretation for X) and allow it to decay 
sufficiently quickly via some suitable modification to this model, the couplings a and/or Kg 
can be freely dialed down without overdosing the universe, thereby alleviating the direct 
detection constraints and allowing light m^. 

In the EFT description for this model, the operator leads to a “Higgs 

oblique” correction to the Higgstrahlung process; i.e., the operator induces a Higgs wave- 
function renormalization 6Zh = which modifies the cross-section by /S.azhl^Zh = 

5Zh [26, 49]. With a sensitivity of ^azh/f^Zh ~ 0.5% (corresponding closely to the CEPC 
and FCC-ee sensitivities in table 1), values of mgja < 2.5 TeV could be ruled out at 95% 
confidence, with 5-a discovery reach up to mgja ~ 1.6 TeV [49]. (The limits on mg/a just 
quoted are slightly weaker than those in figure 1 of ref. [50], wherein a combined fit to all 
projected ILC Higgs coupling measurements is considered.) 

In this model, the electroweak precision observables S and T are generated by one- 
loop running of the operators between fj, = mg and ^ = mw, leading to a mixing of the 
Higgs oblique operator with the operators responsible for giving the S and T parameters. 
This one-loop running has been computed in ref. [50], which considers the UV completion 
above in the limit that ng ^ 0 (and with no assumption on the sizes of eg,b). The 
leading contributions to the S and T parameters are found to be log-enhanced, but linearly 
proportional to the high-scale (// = mg) value of the Wilson coefficient of the Higgs oblique 
operator, and suppressed by a loop factor. 
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To summarize the results of ref. [50] , the electroweak precision constraints on the model 
are marginally weaker than those from the Higgsstrahlung measurement, if the most opti¬ 
mistic Higgsstrahlung precision (with FCC-ee or CEPC) is compared to the most optimistic 
electroweak precision constraints (with FCC-ee-t). Current constraints on S and T exclude 
ms!a < 1.2 TeV at 95% confidence. 

In the discussion above, we have neglected the impact of the neutral x held that is 
still present in the EFT given by eq. (3.3). This objection notwithstanding, provided that 
rriy^ 3> E*^ where E* is the energy scale for some process of interest, one can also integrate 
the X out at one-loop, in which case the Wilson coefficient of the operator will 

shift by a subdominant amount; here, the shift will be smaller by a factor of ~ 
than the leading contribution. Therefore turning on ks as a non-zero weak coupling, and 
integrating the x out at the scale /U = m^, will give a subdominant one-loop shift to czh: and 
higher-loop contributions to S and T. The ratio azh ■ S : T should not change significantly, 
thus roughly preserving the relative strengths of the limits when x is added to the theory. 
The foregoing limits on msja are therefore expected to be fairly accurate. They arise as 
an effect caused by the mediator, and are fairly insensitive to the presence of the portal 
coupling per se. 

In summary, the precision electroweak and Higgsstrahlung constraints on this scenario 
are very similar to those of the model in which the x is simply absent, and this case has 
already received extensive attention in the literature. 

4 Singlet-doublet UV completion 

We devote the bulk of this paper to the singlet-doublet UV completion of the fermionic 
Higgs portal operator. The model consists of the Standard Model augmented by the same 
vector-like Dirac fermion SM-singlet x ^s before, as well as a vector-like Dirac fermion 
5f7(2)-doublet E = {C,N) transforming under the SM gauge group as (1, 2, -|-l/2). The 
coupling of these particles to the SM is taken to be 

£ = £sm + ix^X - ^xXX + iFI^E - MpFE - kFHx - kxH'^F, (4.1) 

where — igW^^'E — ig'B^Y, and E = ^cr“. Without loss of generality we may 

absorb the phase of k into the definition of x and/or E, and so throughout we take k to 
be a non-negative real parameter. Without much modification (it would suffice, e.g., to 
replace m^xx m^XrXR + h-C-, with argm^ / 0) it would also be possible to generate the 
CP-odd Higgs portal operator FffFTx^TsXi "w® do not, however, consider this case further, 
and take all parameters to be real. 

In the parameter region where Mp ^ rriy., the heavy doublet E can be integrated out at 
the scale g, = Mp, and the leading correction to the SM Lagrangian is the fermionic Higgs 
portal operator. We will consider constraints on this UV completion both in the regime 
where Mp ^ rriy., constructing an EFT to analyze the low-energy effects, as well as in the 
more general mass parameter space. For general masses, this model does not necessarily 
provide the UV completion to the fermionic Higgs portal operator as it is defined here 
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and our EFT is not valid, so we instead perform direct one-loop computations of relevant 
observables. 

Before discussing the computations and results in detail, it is worth laying out our 
general expectations. The Yukawa-like coupling n between the F, x H fields in eq. (4.1) 
is a hard breaking of the accidental global SU{2)v custodial symmetry [64] of the SM, 
and as such we expect fairly large corrections to the precision electroweak T parameter. 
Additionally, the mass-splitting (i.e., weak iso-spin breaking) in the F doublet which arises 
from the mixing of the neutral fermions N and x after electroweak symmetry breaking 
(EWSB), leads us to expect that there will additionally be contributions to the electroweak 
S parameter. 

Note that, although we do not consider such a case in this paper, custodial symmetry 
could be restored—or broken in a controlled fashion—by augmenting the field content with 
an additional positively-charged vector-like fermion -0 ~ (1, 1, -|-1) with the same Yukawa- 
like coupling to the H and F fields; the mass-splitting |m^ — m^| then controls the degree to 
which the symmetry is broken (being restored in the degenerate-mass limit). There would 
however be additional experimental handles in this case as the new direct coupling of the 
Higgs doublet to electrically-charged fermions would lead to one-loop corrections to the 
/i —)• 77 rate. 

We also expect that, in the model defined by eq. (4.1), the Higgs oblique correction will 
be generated (by closed N-x loops on the Higgs propagator, in the full-theory picture), so 
that even if we were to tune away the large corrections to T by, e.g., the method indicated 
in the previous paragraph, significant constraints would still remain on this model. To be 
clear, however, the Higgs oblique correction is by no means the only contribution to the 
shift in the Higgsstrahlung cross-section which is induced in this model; indeed, non-zero 
S and T parameters themselves will also generate a shift to azh (see ref. [49]), and we take 
all such shifts into account. 

4.1 Effective field theory approach 

It is instructive to begin our analysis by integrating out the new field content, assuming 
it is heavy, and considering the low-energy effects of the dimension-6 operators that are 
generated. In many popular bases of gauge-invariant dimension-6 operators which extend 
the SM (see, e.g., refs. [65-68]), one can immediately read off the S and T parameters as 
the Wilson coefficients of certain operators. In the so-called “HISZ” basis [67], S is simply 
proportional to the Wilson coefficient of the operator , and T is likewise pro¬ 
portional to the Wilson coefficient of the operator the observable corrections to 

the momentum-dependent part of the Higgs two-point function would arise from the opera¬ 
tors D^H\‘^ and in this basis, but the full azh corrections require significantly 

more work to obtain [49]. 

In the following, we first integrate out the F doublet at the “high-scale” = Mp, 
performing the one-loop matching onto an effective theory with the Standard Model field 
content plus the singlet y. We then repeat the procedure, integrating out the y field at the 
“low-scale” /r = rriy. and matching at one-loop onto a new EFT with the SM field content 
only. We assume the mass hierarchy Mp ^ f ~ where Mp ^ is imposed to 
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Table 2. Operators appearing in the effective theory below the scale Mp (see text for sign and 
normalisation conventions). For operators with a % bilinear, we have performed tree-level matching 
at the high-scale Mp up to and including dimension-10 operators in order to be able to compute 
corrections to fairly high order in rriy^/Mp when later matching at /r = but omit these lengthy 
results here as we used only selected terms from them as necessary (see text). Operators containing 
only SM content appear in the one-loop matching both at the high scale = Mp and at the low scale 
/i = The naming convention for the operators with only n Higgs doublets H and m (covariant) 
derivatives is 0(n,m); other operators are named on an ad-hoc basis. We define Wf^i, = 
and = ig'. Note that the operator 0(2,4) is non-standard in the literature and can be 
eliminated using eq. (4.2). 


Tree-level — operators containing a x bilinear 


Name 

Operator 

Name 

Operator 

O 5 

H^HxX 

Om 

Ixi^X - h.c.) 

Oqb 

i(x^x-h.c.) 

O7A 

lH^H{x\Fx + h.c.) 

O7B 

\DpH\^XX 

O 7 C 

5 * 

{{D^H)^D,H-h.c.)x^P'^X 

O 7 D 

i - h.c.) (x7m^X - h.c.) 

etc. 



One-loop — operators with only SM content 

0(2,2) 

\DpH\^ 

0(2,0) 

\H\ 

2 

0(4,0) 

\H\^ 

0(4,2), A 

Udu\H\r 

0(4,2),B 


0(4,2),C 

\H\ 


0(6,0) 

\H\^ 

0(2,4) 

+ Fc. 

O(WVP) 


0(BB) 


0(BW) 


0(Dm) 

Tr 

[D^,W,p\[DP,W-P] 


0{DB) 


0(w) 


^iB) 


^{WWW) 

Tr 



guarantee that the leading correction upon integrating out the F is the CP-even fermionic 
Higgs portal operator; this is not a necessary assumption of the model itself. 

In performing the matching, we find it convenient to match onto the basis of dimension- 
6 operators with SM fields shown in table 2. This basis of operators is not one of the stan¬ 
dard ones and in particular certain operators (e.g., 0 ( 2 , 4 )) ^re redundant in the sense that 
they are related to other operators by SM classical equations of motion (EOM). The redun¬ 
dant operators can be eliminated by field re-definitions whose effect is exactly equivalent to 
making replacements using the classical EOM (see, e.g., refs. [65, 69]); this is correct even at 
the quantum level if we consistently neglect higher-dimension operator corrections. We note 
for now that the operators in table 2 form a convenient and sufficient set to match onto; we 
will return below to the treatment of the various redundancies in order to transform into 
the operator bases most convenient for calculating the S and T parameters, and crzh- 

Upon integrating out F, we find the operators and Wilson coefficients shown in tables 
2 and 3. The only operators of dimension 4 or less appearing in table 2 - 0(2,0)) 0(4 ,q), 
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Table 3. Wilson coefficients of the operators appearing in table 2, evaluated at /r = Mp, the 
“high scale.” Our sign convention is that the Lagrangian term containing any operator of dimension 
greater than 4 which is shown in table 2 appears multiplied by its corresponding Wilson coefficient 
and the appropriate inverse power of the EFT cutoff-scale Mp\ however, for operators of dimension 4 
or less, the sign appearing in front of the operator follows the SM conventions and any multiplying 
power of Mp is absorbed into the Wilson coefficient (see discussion in text). These expressions 
assume a hierarchy of scales Mp ^ m^, and are correct to tree-level for the Ci and to one-loop for 
the B(j). We present results to fairly high order in m^/Mp, although not all of these terms are 
numerically necessary to obtain our results. 


Tree-level — operators containing a x bilinear 


Coefficient 

Value 

Coefficient 

Value 

C5 


C6A 


cqb 


C-7A 


C7B 


C7C 

-k2 

C7D 

-At2 

etc. 



One-loop — operators with only SM content 
Coefficient Value 

-(KV47r2) Ml [1 + mjMp + + m^/Ml 

-b my Ml + my Ml H-] 

3(2,2) +(KVl6vr2) [1-2 myMp - 4 /M^ - 10 /M| - 15 m^/Ml 

-24:myMl - 32myMl + • • •] 

-^(4,0) —[l + Am^/Mp -b 9m^/M|. + 16m^/M|, -b 25m^/M|, -b • • •] 

B(i,2),A -(K^/247r2) [b + Am^/Mp - ^m^/Ml] 

3(4, 2 ),B -{k‘^/A8tv^) [5 - 2m^/Mp - lb my Ml -|- ] 

3(4,2),c -{b^/A8'k‘^) [l - 22 myMp - l22,myMl -|-] 

3(^,0) -{k^IV2'k^) [l + l2myMp + A8ml/Ml ■ ■ ■] 

3(2,4) +(K^/487r^) [l - myMp - 3m‘^/Ml - 20myMl - AO my Ml H- ] 

3(ww) -{b?/IAAtx'^) [l - AmyMp + 3myMl -b 10m^/M|, + OSm^/Ml H-] 

3(bb) -{k^/IAAtt'^) [l - AmyMp + 2>myMl -b 10m^/M|. + 68H-] 

3(bw) -(K^/727r^) [l - Am^/Mp + 'Amy Ml + 10 my Ml + 08ml/Ml ^-] 

3(dw) +(l/2407r^) 

3(db) +(l/2407r^) 

3(w) +{K^/72Tr^) [7-4 m^/Mp - 12 my Ml -|- ] 

3(B) /7‘2.'k‘^) \J - Am^/Mp - l2myMl a -] 

3(WWW) Operator generated, but coefficient not computed. 

Irrelevant to the observables of interest. 
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and 0 ( 2 , 2 ) ^ already appear in the SM. Our convention for these operators is such that the 
Wilson coefficients - 6 ( 2 , 0 )) 6 ( 4 ^ 0 )) ^{ 2 , 2 ) simply add positively to the coefficients already 

present for these operators in the SM: i.e., £ D (1 + 6 ( 2 ^ 2 )) 0 ( 2 , 2 ) + + 6 ( 2 ,o)) 0 ( 2 ,o) — (A + 

6(4,0)) 0(4^0)- the other hand, all operators with dimension greater than 4 are taken to 
appear in the Lagrangian as £ D Cj Mp~ * Ot where Cj is the Wilson coefficient and di is 
the dimension of the operator Oi. 

We worked to tree-level for the Wilson coefficients for any operator containing a x- 
field bilinear, and to one-loop in obtaining the Wilson coefficients for operators containing 
only SM field content. The rationale is that whenever the x later appears in any diagram 
computed in a consistent one-loop matching to obtain a Wilson coefficient of a dimension -6 
operator containing only SM content, it will appear as a closed x loop, so only the tree-level 
X-SM couplings are required. 

The tree-level matching was performed by inverting the classical equation of motion 
for the heavy doublet to solve for F in terms of H and y, and expanding in powers of 
D^/Mp (some lower-dimensional results were checked diagrammatically). The one-loop 
matching was performed diagrammatically utilizing dimensional regularization and the MS 
renormalisation scheme. Calculations were done in the broken electroweak phase, and on 
the EFT-side of the one-loop matching computations we consistently included diagrams 
with one or more tree-level X'SM coupling(s) and a single closed x loop. As a check, some 
one-loop computations were performed both in the non-mass-diagonal basis of new neutral 
fermions, N and X) and in the mass-diagonalised basis. ^ In the latter case, expansion of 
the result in powers of the Higgs vev v allows the comparison of these two methods and 
the extraction of the Wilson coefficients. Finally, we have checked that, to the order we 
work, the EFT and full theory one-loop computations reproduce the same IR-divergent^ 
logarithms (~ log at fi = Mp) as is of course required. 

The EFT power-counting we employed is such that we kept only leading terms in v/Mp, 
but retained terms at higher order in m-^^/Mp, although this was not always numerically 
necessary to obtain the results we present below. In order to compute these higher-order 
terms in the mass-ratio m^jMp in our one-loop matching, we have systematically carried 
out the tree-level matching up to and including dimension -10 operators so as to capture 
the operators with higher numbers of derivatives acting on the x held since, in dimensional 
regularization and MS, such higher-derivative operators contribute positive powers of 
to the S'-matrix element when appearing in conjunction with a closed x-loop. This latter 
point is perhaps opaque, and is illuminated by a simple concrete example: the term in 
Bfypw) suppressed by m^/Mp relative to the leading term can be found from, e.g., the 
v‘^{p^/Mp){rn^/Mp) part of the ZZ two-point function at one-loop (along with other VV 

^Mixing-angle effects are necessarily small in the range of EFT validity, being suppressed by at least 
kv/Mf for Mf ^ m^,KV (see eq. (A.7)). It Is well-known (see, e.g., §3.3 (pp. 37-38) of ref. [70], and 
references therein) that physical effects which In the full theory would be ascribed to these small mixing 
angles are instead captured in the EFT description through the inclusion of higher-order operators, whose 
Wilson coefficients have their genesis (in part) in the UV-theory mixing angles. Readers unfamiliar with 
this point may wish to consult Appendix B for further comments. 

^In the sense that they diverge if —>■ 0. 
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two-point functions so as to break the degeneracy between and B(^bb))- 

Simple power-counting then indicates that, owing to the suppression, this requires in 
turn knowledge of the subset of the operators at dimension-10 in the tree-level matching 
which contain two Z boson fields, two Higgs fields, one x held, one y field, two derivatives 
acting on the Z fields, and-to be dimensionally correct-one derivative acting on either of 
the X or X fields. When the x loop is closed and the Higgs field replaced by its vev, such 
operators contribute to the ZZ two-point function in the EFT with a term proportional to 
{v‘^p^/Mp) f d'^kk'^/{k'^ — m^) ~ {v'^p^/Mp)m‘^, as desired. 

Having obtained the mixed tree-level/one-loop matching conditions, one should by 
all rights run the theory down from the high-scale p = Mp to the “low-scale” p = 
integrate the x out at this scale, and then run the theory to the scale required to consider 
any process of interest (i.e., the matching conditions supply only the initial conditions 
for the renormalisation group equation (RGE) running). However, we will choose to not 
run the operators using the RGE. The errors one makes in ignoring the running between 
pL = Mp and p = are proportional to log provided we do not take too 

large a hierarchy, this factor is not necessarily very large. Furthermore, for the electroweak 
precision observables S and T, we expect the effect of the running is either higher-loop or 
power-suppressed in m^/Mp and therefore a subdominant modification to the coefficients 
generated at the high-scale. This is in contrast to, e.g., the case considered in ref. [50], 
where the operators responsible for S and T are not generated at the high-scale, and arise 
purely from running: in that case, it is important to find the contributions from running 
to the corresponding Wilson coefficients because S and T constraints are so strong. Our 
neglect of the RGE in this model is also numerically justified, at least for the S and T 
parameters, as we will show in comparing EFT and full-theory one-loop results. 

We thus move on to the matching at the low-scale fj, = rriy. where we integrate the x 
out, working consistently correct to one-loop in the matching. The operators Oi for i G 
{5, 6H, 6B, 7A, etc.} are thereby removed from the Lagrangian and the Wilson coefficients 
Bi for the dimension-6 operators with only SM field content receive new contributions in 
the EFT valid below /x = rriy’. We call these shifted Wilson coefficients in the low-scale 
EFT Ci] they are shown in table 4 assuming no running or operator mixing between the 
scales p = Mp and p = niy-. 

While the basis of operators shown in table 2 was convenient for the matching com¬ 
putation, it still contains the various redundancies we mentioned earlier. It also does not 
happen to be the most convenient for the computation of the S and T parameters, and 
(Jzh- The S and T parameters are each simply related to the Wilson coefficient of a single 
operator when these redundancies are removed. In addition, ref. [49] contains a computa¬ 
tion of azh in terms of a different basis of operators, and it useful to transform to this basis 
to enable us to perform a cross-check on our own independent computation of azh- 

As d( 2 , 4 ) is non-standard in the literature, we eliminate it via the identity 

^*(2,4) = “ (^{ww) — 20(5) — 20(vj/), (4.2) 

plus irrelevant total-derivative terms. We are also free to make use of the SM EOM for 
H to trade out for corrections to the coefficients of the operators |//p,|i/|^,|i^|® 
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and the SM Yukawa Lagrangian in addition to two new types of dimension-6 

operators —H x and a number of four-fermion contact operators proportional 

to two powers of the Yukawa couplings—as well as dimension-8 corrections, the latter of 
which we consistently neglect. The corrections to the and terms are 

absorbed into unobservable shifts of the SM parameters and Tj (the SM Yukawa 

coupling matrices), but the remaining additional terms here lead in principle to observable 
effects. 

Our last step is to rescale the H field to obtain a canonically normalized covariant 
derivative term (i.e., we absorb the coefficient of (1 + 0(2,2))) iiito the H field) 

and re-define the SM couplings [;U^, A, and Tj] to absorb both the shifts proportional to 
C( 2 , 2 )) 3-8 well as the shifts to the SM parameters arising from 0(4 0 ) ^( 2 , 0 )- We note 

that in the coefficients of all dimension-6 operators, the rescaling of the H field can simply 
be ignored as the leading coefficients are already linear in the O^: any correction resulting 
from the H rescaling is therefore proportional to 0(2,2) ^ Oj and thus formally counted as a 
two-loop effect, since a single Oj insertion is formally counted as a one-loop correction. In 
other words, the coefficients 0 ( 2 , 2 )) ^*( 4 , 0 ) and 0 ( 2 , 0 ) can simply be ignored in an analysis 
correct to one-loop (this assumes a degree of tuning in the bare value of such that we 
still obtain the correct measured Higgs mass). 


4.1.1 Precision electroweak observables 

Having eliminated 0(2,4) using eq. (4.2), the S and T parameters can be read off from the 
new coefficients of the operators 0(5^( and 0(4 2 ),S) respectively.^ These operators modify 
the gauge boson vacuum polarization amplitudes Yi^y{q^) = + ••• • Correct 

to linear order in the Wilson coefficients Oj, we find that the new physics contributions to 
S and T are:® 


4 p2 „2 


Ote 

.2 „2 


TlzziO) 


Ote 

Acwsw 


nixo) - 

civsiv 


CIV'S TV 


M| 


n'2,(o) - ng(o) 


n^y3_B(0) 


[‘^(StV) - 20(2,4)] ) 


M| 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


®Again, we do not run the operator coefficients between and the Z-pole. 

®The operators 0(db),{dw) contribute only terms ~ to the gauge-boson two-point functions and 
hence do not contribute to the S, and T parameters (at least if S is defined per eq. (4.4)); instead, these 
operators give rise to the Y and W operators of ref. [71]. However, we find the corresponding limits are 
not competitive with those from S and T and do not pursue this further. 
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Table 4. Wilson coefficients of the operators appearing in table 2, evaluated at ^ = rriy^, the “low 
scale,” after the x has been integrated out. These expressions are correct under the assumption of 
no operator running or mixing between the scales ^ = Mp and ^ = m^. Our sign conventions are 
as detailed in table 3. These expressions assume a hierarchy of scales Mp ^ m^. 

Coefficient Value 


C( 2 ,o) — (K^/47r^) \_Mp + Mp + m^] ( this appears to be exact ) 

^2,2) +(KVl6vr2)[l - 2mjMp - 

-10 ml /Ml -12 ml + ''' ] 

C'(4,0) [l + im^/Mp + 8ml/Ml + 12ml/Ml + lOml/Mj,-\ -] 

C'(4,2),A -(K^/247r2) [7 + 8m^/Mp - 8ml/M], -] 

<^(4,2),^ -(K^/487r2) [5 - 2m^/Mp - 15H-] 

C'(4,2),c -(K^/487r^) [l - 22 m^/Mp - 111 ml/Mj. H-] 

^(6,0) -(K®/127r^) [l + lAm^/Mp + 51 H-] 

C{ww) -(K^/1447r2) [l - Am^/Mp + 3m^/M|. + 10ml /+ 41 H-] 

C(BB) -(KVl447r2) [1 - AmjMp + + 10 my + Alml/Ml + • • •] 

C{BW) -(K^/727r^) [l - Am^/Mp + 8myMl + 10+ 41 ml/Mj, H-] 

^(2,4) +(K^/487r^) [l — m^/Mp — 3m^/M|. — 14m^/M|, — 25m^/M|, + • • •] 

C{db) +(l/2407r^) 

C(DW) +(l/2407r^) 

C'(W) +(K^/727r^) [7-4m^/M^- 12 myMj, A -] 

C(B) +{K^/72n'^) [7 - Am^/Mp - 12myMl A -] 

C{WWW) Operator generated, but coefficient not computed. 

Irrelevant to the observables of interest. 


and 


T = 


P- 1 
CXe 




^eM^ L 


n^y+p/-(0) — Ylyp3y/z{0) 

71w+w- (0) — c^n2^(o) 


— —c 
2a^Ml 


(4.7) 

(4.8) 

(4.9) 
(4.10) 


where in arriving at Eqs. (4.4) and (4.9), we have made use of liz^O) = n.y.y(0) = 0. The 
U parameter is identically zero at dimension 6; the hrst non-zero contribution occurs at 
dimension 8. 
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Substituting the explicit forms of the Wilson coefficients, we have 


2k^ v‘^ 

97 


1 - T 


7 3 rrii mt 17 rnt 

A. _^ Q A _^ 

2 Ml Ml 


4M, 


2 Ml 


+ 


(4.11) 




Oe 247r Ml 


2 

5 Mp 


'Ml 


where = — (4-12) 

47r 


?7 = 0 + (dim-8) (4-13) 

The parametric enhancement (or suppression, depending on the relative sizes of and 
Ue) of T compared to S can be understood as follows. For Mp S> m^, ^y[/+^y-(0) — 
cl^UzziO) ~ 7(7n^)(7n^/M|.), while 1122 ( 0 ) ~ {g^ + (7)^)(7n^/M|.). This is deter¬ 
mined by the mass dimension of the quantity concerned and the fact that, in our model, 
a diagrammatic new-physics contribution proportional to the Higgs vev v must always ac¬ 
tually be proportional to kv owing to the structure of the H-F-x coupling. Therefore, 
T ~ (77/M^)(7/ae)(77/M|) ~ {ai^lae){K'^v‘^jMl) while S ~ (c^s^/ae)n'22(0) ~ 
k^v'^/MI since ~ 5^ + parametric dependence on k can also be read 

off from the form of the dimension-6 operators giving rise to S and T. Upon EWSB, the S 
parameter arises from O(biv) = H. As a result, diagrammatic contributions to 

S come with two powers of the Higgs vev v and hence two powers of nv, since the two deriva¬ 
tives in the field strengths are removed in the definition of (S' as a momentum derivative of 
the two-point function. The T parameter arises, upon EWSB, from 0(4, 2 ),B = 
the diagrammatic contributions to T thus come with four powers of the Higgs vev v and 
hence four powers of kv (two of the Higgs vev factors are later cancelled by the M^ in the 
definition of T). Therefore, T ^ k'^ and S ^ k?. 

Note that an indirect effect of the dimension-6 operators may be to imply shifts to the 
Lagrangian coupling constants (e.g., the gauge couplings) away from their SM reference 
values, as we will discuss when we turn to the Higgsstrahlung process in the following 
section. However, since the leading results for S and T are already proportional to C* and 
hence formally one-loop, any coupling shift can be dropped as a higher-order correction 
((S' = T = 0 in the SM at tree-level). Moreover, the coupling constant shifts will imply a 
shift in the numerical values of the one-loop computations of the SM contributions to the 
S and T parameters, but these are again formally at least two-loop effects. Since we work 
to one-loop accuracy, we neglect these contributions and the value of appearing in T is 
taken to be its SM value (we take it at the Z pole). 

In figure 2, we show the 95% confidence exclusion regions in the mass parameter space 
from measurements of the electroweak precision variables (5, T), for fixed representative 
values of k. These results are computed in the EFT, and are shown for both current limits 
(LEP+SLD) and for the proposed CEPC, ILC and FCC-ee colliders with the sensitivities 
given in table 1. 

In the region of parameter space under consideration in figure 2, the limits are driven 
almost exclusively by the size of the T parameter, owing to its parametric enhancement 
~ Ui^jae over the S parameter. In the regime where the EFT is valid, < \Mp, the 
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Figure 2. Two-parameter 95% confidence exclusion regions (—21og[£/£o] ^ 5.99) from measurement 
of the precision electroweak variables (5, T), using results computed in the EFT. We present these as 
boundaries in the allowed mass parameter space for fixed representative values of k, as annotated on each 
line; the unshaded region to the lower-left of these lines is excluded for the given value of k. To compare EFT 
results with the full-theory calculation, see figure 4. The light shaded region, Mp/4 < < Mp/2, denotes 

the region where the EFT begins to break down: the error in the EFT result for T compared to the v'^ /Mp 
piece of the full result (i.e., the “dimension-6 part of” the full result) is ~ 20% at ~ Mf/4, reaches 
~ 50% at ~ Mp/Z, and becomes > 100% before m-^ ~ Mf/2. In the dark-shaded region, m-^ > Mf/2, 
the results have consequently been masked as they are invalid. The various line styles correspond to the 
current constraints and various projected constraints on [S, T) for the proposed CEPC collider (left plot) 
and the proposed ILC and FCC-ee colliders (right plot); see table 1. 


exclusion reach is also largely insensitive to the value of as one might expect. Already 
with current (LEP+SLD) constraints on S', T, masses Mp below about 675 GeV can be 
ruled out for k ~ 1; due to the dependence of T, this lower limit increases to around 
2.9 TeV for a coupling k ~ 2, but there is essentially no limit in the regime where the EFT 
is valid if k 0.3. 

With even the “baseline” proposal for the CEPC, these lower limits increase by a factor 
of approximately 1.75 owing to the factor ~ 3.5 times stronger limits on S and T in this 
scenario as compared to the current bounds [43]. For the CEPC scenario with improved 
measurements of P^ and sin^ 6 only, the lower limits on Mp increase by a factor of ~ 2 
compared to the present limits, rising to ~ 1.5 TeV (6 TeV) for k ~ 1.0 (2.0), in the limit 
where Mp. This results from the tightening in the limits on T (and S) by factor 

of ~ 4 compared to the present constraints [43]. For the best-case scenario for the CEPC, 
with improved measurements of P^, sin^ 0, and rrit, the lower limits on Mp increase by a 
factor of ~ 2.7 compared to the present limits (i.e., about 30% over the previous scenario), 
rising to ~ 1.8 TeV (7.7 TeV) for k ~ 1.0 (2.0), in the limit where <C Mp. This results 
from the tightening in the limit on T by factor of ~ 7, and the limit in S' by a factor of 
~ 5 compared to the present limits (i.e., the mt measurement roughly halves again the 
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measurement uncertainty on T compared to the previous scenario, and also improves the 
S constraint) [43]. 

The projected sensitivities of the ILC and FCC-ee show improvements in the lower 
limits on Mp which are broadly similar to the various CEPC scenarios. The ILC projection 
and CECP “baseline” projection are very similar, as are the FCC-ee-Z limits and the CEPC 
limits with improved measurements of P^ and sin^ 9 only. Finally, the FCC-ee-t lower-limits 
on Mp (for k ~ 1.0 and k ~ 2.0 they are 2.1 TeV and 8.5 TeV, respectively) are slightly 
stronger than the best-case CEPC scenario, but are nevertheless still broadly similar. 


4.1.2 Higgsstrahlung 

The computation of shifts to the Higgsstrahlung cross-section azh is more involved than 
that required to obtain the S and T parameters. There are a large number of operators 
which directly contribute to a shift in the cross-section, both by way of contributions to 
wavefunction renormalization of the h, Z fields in the hZ^Z^ coupling, and the introduction 
of additional diagrams. Furthermore, azh is nonzero at tree-level in the SM and so the 
one-loop new-physics shifts to the relationships between SM input parameters and SM 
Lagrangian parameter values must be accounted for here to maintain a result consistently 
correct to one-loop order: the standard set of SM input parameters for high-precision work 
are {Gp,mz,ae) and the dimension-6 operators impact the processes (e.g., muon decay at 
= 0) used to relate the numerical values of these parameters to the values of Lagrangian 
parameters, leading to shifts in the Lagrangian coupling constants away from their SM 
reference values [72]. 

The requisite computation of azh has actually recently appeared in the literature in 
ref. [49]. We have nevertheless independently repeated the computation of azh hi the EFT 
as a cross-check; we find complete agreement with their results. In order for us to make 
reference to the results of ref. [49], we transform to the same basis of operators used there. 
In addition to the elimination of 0(2,4) described above, four further sets of manipulations 
on our basis of operators are required: (a) we re-write the operator 

0(4,2),B (4.14) 


(b) we rewrite the operators and 0(£)iy), up to total derivative terms, as 


^(DB) 

^(DW) 


iff 


/^2 


{d^B,p){d^^B''P) = and 


Tr 


[o^,LF,p][o^ly"^] 


= 2Tr 


[Dp,W^^P][D,,W^ 


-F 4Tr 


W ^IF PW P 

** p U ** p 


(4.15) 

^ p 

p J ’ 

(4.16) 


and utilise the SM EOM for B^i, and to rewrite the first term on the RHS of each of 
these relationships; and (c) we make two manipulations to remove 0(4 2 ),c- 




(4.17) 
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followed by the substitution of the SM EOM to rewrite D^H. Finally, (d) we rewrite the 
operators 0 ( 5 ) and 0(n/)) up to total-derivative terms, as 


0 ( 5 ) 

0(vy) 


2^{bb) + 2^(bw) - 

2^{ww) + 2^iBW) 


4 

5 
' 4 


hUd^h 


, and 






(4.18) 

(4.19) 


and then utilise the SM EOM for and W^i, to re-write the last term on the RHS in 
each line, and the SM EOM for H to re-write some factors of D^H which appear upon 
doing so. This is followed by a final re-definition of the SM Higgs-Lagrangian parameters, 
to absorb some unobservable shifts to the SM parameters which occur during this process. 
When the dust settles, we are left with 

i ^ 

where the dimension-6 operators and Wilson coefficients C,-, correct to linear order in 

J 

the original C*)...), are listed in table 5, and Csm takes the same form as the SM Lagrangian, 
except that all the parameters are now understood to be defined so as to absorb any of 
the unobservable shifts which occurred as a result of the manipulations just described.^ As 
far as the e"’“e“ —?• Zh process is concerned, all of the operators in table 5 listed above the 
horizontal line contribute, either through a shift in the SM input parameters or through 
the addition of a new term in the amplitude itself.® 

The results of ref. [49] can then be used to read off the value of azh] iu table 6, we supply 
an explicit dictionary to transform our Wilson coefficients into the notational conventions 
of ref. [49]. In all our results, we utilize the {Gp, Mz^o-e) input-parameter set. 

Apart from the independent re-computation of azh beginning from the basis of opera¬ 
tors and Wilson coefficients in table 5 to cross-check against the results of ref. [49], we also 
verified the correctness of the fairly complicated manipulations which were used to remove 
the operator 0(5iy) (see eq. (4.16) and the text following) in arriving at eq. (4.20). To do 
so, we augmented the SM Lagrangian with 0(5VK) alone and considered its effects on azh 
directly, rather than eliminating this operator by EOM. 

Furthermore, since the ILC measurement projections assume a polarized electron beam 
with (Pg_,Pg+) = (—0.8,-1-0.3) [44], in the course of our independent computation, we 
also computed the polarised-beam cross-section. In the notation of ref. [49], the polarized 

^These unobservable shifts can consistently be ignored; they are not the same as the physically relevant 
shifts which need to be accounted for in relating the input parameters to the values of the Lagrangian 
parameters. 

®Note that C(wvc) = C'(sb) = |C'{sw) and C(w) = C'(s) in our model, so the Wilson coefficients of the 
operators H and H"' H remain equal, and equal to one-half of the Wilson coefficient 

of the operator H‘' H. This ultimately follows from the fact that there are no new Higgs-charged- 

charged vertices in the full theory picture, and guarantees that there is no leading-order (one-loop) correction 
to the h —>■ 77 decay rate. This would not be true if we had introduced a charged partner to the x to 
maintain custodial symmetry, as we discussed earlier. 
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Table 5. Operators and corresponding Wilson coefficients Cj appearing in the effective 
Lagrangian which we utilise to compute azh [see eq. (4.20)]. Only the operators listed above the 
horizontal line contribute a shift at leading order to e“'"e“ —>■ Zh, either through a shift in the 
relationships between the SM input parameters and the Lagrangian parameters, or through the 
addition of a new term in the amplitude itself. A, p. and the Yukawa matrices are understood to 
take their unobservably shifted values. We define Wfj,„ = = ig'— Dfi 

and cr“. Superscript p,q are generation indices; i,j are S'[/(2)-fundamental 

indices; and a is an S'{7(2)-adjoint index—all are summed over if repeated. 


Operator Wilson coefficient Cj 



2^2 

“^7(4,2),C + 5^7(4,2),B + C’(4,2).yl “ |S'^C7(W) 



+ ^9^^{DW) + fff^C(2,4) 


~C^(4,2),B — i(S'0^C'(S) + {9')^C(jjb) + (ffO^C'(24) 


C{WW) + k^(W) - 20(24) 


C(SB) + k^iB) ~ 20(24) 


C(BW) + 1 (C(W) + C'(B)) “ 40(24) 


“(flO'^C'pB) + 3(ffO^ ~ 2C’(24)) 


-^9^C(^B)W) + (C’(iv) “ 20(24)) 



\H? 


+C(6,0) + -^^80(24) -1- 20(44)_c — ‘ig'^Ci^BW) 



+5^C'(W) “ 2p^O(24) j 

Tr 


(d(WWW) + 40(£)iy) 

{H^H) 

' {uree)w + {Qirdd)Hi' 
+ {Qiruu)ejk{H^f 

+ h.c. jC(^ 4 ^ 2 ),c ~ ^g^C(DW) + \9‘^C!(W) + 5 (8A — g'^) 0(24) 


-gS'"‘C'(mv) + (C(iv) “ 20(24)) 


- (sO'^c’pB) 


— 49 ^(DW) 

{Qlra-Ql) 

-29 ^(DW) 


{LV,e){eVjL) + {QTad){dVd^Q) 

+{QVuu){uVu^Q) 

+{{lV,e){dVa^Q)+\i.c.) 

-((LWee)e,fe(Q'=r„M)+h.c.) 

-((gWrfd)e,fe(Q'=r„u)+h.c.) 


20(24) 


cross-section is obtained from the unpolarized cross-section by making the following simple 
replacements in the quantities defined either in their eq. (2.3) or in their table 2: multiply 
every appearance of g\ in their F^u and as well as every appearance of gp in their 
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Table 6. A dictionary to convert our Wilson coflicient results from table 4 to the Wilson coefficients 
defined in ref. [49]. 

Wilson Coefficient in ref. [49] Value per our table 4 

“I + C'(W) “ 4C(2,4)) 

“I + C'(B) “ 4(7(2,4)) 

“I (^^(BW) + C'(TT) + C'(B) - 8(7(2,4)) 

“(^(4,2),C + 5 (^( 4 ,2),B + (^(4,2),A 

+ ifi''^(^(BVK) “ (^(W) “ 2(7(2,4)) 

- 5 ^( 4 , 2 ),B + ((^(B) “ 2(7(2,4)) 

-l9*C(^DW) + (C'(vy) - 2(7(2,4)) 

“iff (-"(BVK) 

1{9')‘^C{db) - li9')‘^ ((^(B) - 2(7(2,4)) 

(fl'0‘^(^(DB) - 1 ( 5 ')^ ((^(B) “ 2(7(2,4)) 

^ 3 , 4 , by (1 — Pg-)(1 + Pe+). Similarly, multiply every appearance of g]^ in their F^m and 
Fi, 2 , as well as every appearance of gn in their i^ 3 , 5 , by (1 + Pg-)(1 — Pg+). No changes are 
needed to the factors appearing explicitly in their eq. (3.9). 

In hgure 3 we show the regions in the mass parameter space (m^, Mp) which, for 
the given hxed value of k, would yield a value of azh in conflict with the projected 95% 
conhdence limits as given in table 1. We can see that the best sensitivity comes from the 
CEPC and FCC-ee. At the CEPC, the lower bounds on Mp range from 590 GeV to 7 TeV 
for K ~ 1.0 — 4.0® in the limit where <C Mp. These lower limits rise to approximately 
660 GeV to 10 TeV for ~ Mp/A, which is roughly where we begin questioning the 
validity of the EFT. (The limits on Mp in the region where the EFT results are valid 
are more sensitive to the value of than the electroweak precision limits.) The FGG-ee 
constraints are again a little stronger than GEPG constraints. Lower limits on Mp in the 
limit where rriy’ <C Mp range from approximately 630 GeV to 7.8 TeV for k ~ 1.0 — 4.0. 

The constraints obtained at the ILG are somewhat weaker in all three of the considered 
scenarios than either the GEPG or FGG-ee constraints. Note that care was taken here in 
combining limits where runs are at different energies as {Aazh)l^^ is energy-dependent 
(see comment in caption on table 1). The most optimistic ILG scenario with several ab“^ of 
data yields lower limits on Mp for <C Mp that are approximately 490 GeV and 5.9 TeV, 
respectively, for k ~ 1.0 and 4.0. 

In none of these cases are the limits from the precision Higgsstrahlung measurement 
competitive with the electroweak precision programs at these future colliders in imposing 
constraints on this specihc model; nevertheless, these results do demonstrate that the azh 
measurement would provide a strong complimentary constraint on closely allied models 
where the T parameter is dialed away, as we discussed previously. 

= 4.0 is a fairly large coupling; we display these results only with the explicit point of indicating that 
a very large coupling is required to probe this region of parameter space. 


Cww 

cbb 

CWB 

CH 


Cp 
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Figure 3. One-parameter 95% confidence exclusion regions (—2 log[>C/£o] ^ 3.84) from precision 
measurements of azh, using results computed in the EFT. We present these as boundaries in the 
allowed mass parameter space for fixed representative values of k, as annotated on each line; the 
unshaded region to the lower-left of these lines is excluded for the given value of n. Absent the full 
loop computation, it is not possible to quote an error on this EFT-based result, but based on the 
comparison of the EFT and full-theory computations for the EWPO results, it is probable that the 
EFT results here are questionable in the light shaded region, Mp/i < < Mp/2, and are almost 

certainty invalid in the dark-shaded region, rriy, > Mp/2, where the results have consequently been 
masked. See table 1 for the assumed sensitivities of the experiments. 


Although a large number of operators contribute to azh, a partial and heuristic un¬ 
derstanding of the generic difference in the strength of these limits can be obtained by 
examining a subset of the operators relevant for the generation of T and Aazh- Sup¬ 
pose L = £sm + — h.c.)'^ + . It can then be easily shown that 

T = (au^)/(aeA^) ~ 10^o(u^/A^), and a little more work shows that Aazhl<^Zh ~ 
— (u^/A^)(6 -|- 0.83a);^*^ note that the two operators contribute almost equally to the Hig- 
gstrahlung cross-section. The anticipated one-parameter 95% confidence measurement un¬ 
certainties on T (restricted to S = U = 0) are 2.0 x 10“^, 1.5 x 10“^, and 8.9 x 10“^ for 
the CEPC baseline, “Improved r2,sin^0”, and “Improved Fz, sin^ 0, mt^ scenarios, respec¬ 
tively.^^ The resulting 95% confidence lower bounds on A/-^/|a[ are approximately 20 TeV, 
23 TeV and 29 TeV, respectively. On the other hand, the 95% confidence measurement 
uncertainty on the Higgsstrahlung cross-section at CEPC is projected to be ~ 1% (corre¬ 
sponding to the 68% confidence projection Aazh!(^Zh — 0-5% in table 1), which yields the 
limit A/-|- 0.83a| > 2.5 TeV. It is clear that the latter bounds are significantly weaker 
than the EWPO constraints for roughly equally-sized Wilson coefficients, which is the sce- 

^°The ‘0.83’ here is the numerical value of a fairly complicated function of the gauge couplings, assuming 
g = 0.648,5' = 0.358. 

^^These are obtained from the values in table 1 as = 1.96 ut -^/l — p|-r- 
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Figure 4. As for figure 2, two-parameter 95% confidence exclusion regions (—21og[£/£o] ^ 5.99) 
from measurement of the S and T parameters, except these results are computed using the full one- 
loop results described in section 4.2 and hence are valid for the full range of masses shown. Note 
that and Mp are the input mass-parameters, and are not the physical masses of the neutral 
fermions; see Appendix A. 

nario one would expect when these operators are generated at the same loop order (provided 
of course that there is no symmetry—e.g., custodial—which prevents this). Parametrically, 
it is clear that the relative strength of the limits can be traced directly to the enhancement 
of T by a factor of l/ag. 

4.2 Loop computation of EWPO 

By design, the EFT we have constructed is valid only in the region of parameter space 
where Mp ^ and it is in this region of parameter space where this model supplies a 
UV completion of the fermionic Higgs portal operator H"' HxX- could of course also 
construct other EFTs valid in the region m-^,Mp 3> ~ fnw^ but for similar particle 

masses, ~ Mp, or for the opposite mass hierarchy, my. 3> Mp. Given the amount 
of work required to obtain the results already presented, it appears to be more efficient 
to simply perform the full one-loop computations of S', T; this computation is detailed in 
Appendix A and was performed in the mass-eigenstate basis of neutral fermions, accounting 
for mixing angle effects. The full loop computation of azh is beyond the scope of this paper, 
as we have already established from the EFT results that the Higgsstrahlung bounds are 
expected to be somewhat weaker than those from EWPO. 

The exclusion regions due to precision electroweak constraints arising from the full 
loop computation of S, T are shown in figure 4. In the limit where <C Mp it is clear 
that the results agree well with the EFT computation (c.f., figure 2), the small differences 
being ascribable to our neglect of dimension-8-and-higher operators, as well as the running 
between scales, in the EFT. Note that if the full result for S, T is expanded out to 0{v‘^/Mp) 
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to find the contribution to the result ascribable to dimension-6 operators, we find all the 
same analytic (in the masses) pieces of the result as in the EFT computation, but naturally 
also find non-analytic factors ~ log[m^/Mj,] which are absent in our EFT result as we 
neglected the running of the Wilson coefficients between scales ^ = Mp and ^ 
however, for the contributions to S and T, these are suppressed by at least the third power 
of m^/Mp and are thus negligibly small whenever the EFT which we constructed is valid. 

As we have already discussed the results in the region <C Mj? in connection with 
the EFT results, we focus here on the parameter space not covered by the EFT. Of course 
in this region, the theory does not provide a UV completion of the fermionic Higgs portal 
operator H^HxX: although in the limit ^ Mp^ it would provide a completion of the 
fermionic-doublet Higgs portal operator H^HFF. 

The most obvious point is that a significantly larger region of the rriy- parameter space 
at small Mp is ruled out than vice versa. For example, the fully improved CEPC results 
indicate that if Mp = 100 GeV, then niy, < 33 TeV can be constrained for k ~ 2.0, 
compared to Mp < 7.7 TeV being constrained if = 100 GeV. This pattern is generic 
for all of the experiments; it traces its origin to the fact that for <C Mp, S and T 
are both positive and, roughly speaking, IS"! ~ 0.1|T| near the exclusion limit, whereas 
for Mp m^, T is positive and S is negative, with IS"! ~ |T| near the exclusion limit 
(indeed, IS*! > |T| is possible here). The larger deviation from {S,T) = (0,0) in this region 
leads to the stronger limits. There is also a stronger dependence of the boundary of the 
exclusion region on Mp sX large than vice versa. Both of these behaviours are due to 
the fact that, parametrically, S ~ AM/Mp where AM is the mass-splitting in the doublet: 
AM = = v{k‘^/ 2)[v/{Mp — m^)]. This parametric dependence arises 

because S is zero in the weak-isospin-symmetric limit of the theory, and because it is always 
F fermions which couple directly to W and Z. This implies that the leading contributions 
are S ~ K?v‘^/Mp for <C Mp, but S ~ —K^v‘^/{Mpm^) for Mp <C clearly then, at 
light Mp and heavy m^, IS"! will be larger than if the mass hierarchy were reversed. 

The slightly different shapes of the ILG and FGG-ee exclusion regions when rriy. > Mp, 
as seen in the right plot of figure 4, arise from the different alignments of the 68% coverage 
likelihood contours in the {S, T) plane [43] for the ILG compared to those for FGG-ee, which 
in turn is due to the different projected sensitivites to mw, sin^ and T^ at each collider 
(see also ref. [73]). With reference to figure 1 in ref. [43], it is clear that if S', T > 0 and 
jSj <C |T|, as is the case for <C Mp, the limits from FGG-ee-Z (TLEP-Z) are more 
constraining than those for the ILG; whereas if T > 0, and S ~ —T, as is the case for for 
Mp rriy-, the ILG and FGG-ee-Z limits will be about equally constraining, with the ILG 
actually marginally better, as is visible in the right panel of our figure 4 for Mp < 150 GeV. 


5 Conclusions 

In this paper we examined two possible models which UV-complete the GP-even fermionic 
Higgs portal operator and investigated how future precision electroweak and pre¬ 

cision Higgsstrahlung measurements at the proposed ILG, FGG-ee and GEPG high-energy 
e'’'e“ colliders could constrain these models. In the first model, the “scalar completion,” 
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a scalar singlet acts as a mediator between vector-like Dirac fermion x and the SM Higgs 
field H, with a Yukawa coupling ksSxx to the x and a coupling amsS\H\‘^ to the Higgs 
doublet. In the second model, the “fermionic completion,” a vector-like SU{2) doublet of 
Dirac fermions couples to the x and H field with a Yukawa-like interaction kFHx + h.c.. 

For the scalar completion, we estimated that the effect of the new x particle will 
be subdominant (i.e., loop-suppressed) for the computation of the electroweak precision 
variables S and T and the Higgs couplings to the SM particles, compared to the dominant 
effect of the singlet mediator itself. The limits which can be placed on the SM augmented 
with such a singlet scalar mediator have already been extensively studied in the literature 
(e.g., ref. [50]): a precision Higgsstrahlung cross-section measurement {Aazh)l(^Zh ~ 0.5% 
can place a 95% confidence upper limit on rnsja around 2.5 TeV, while the precision 
electroweak limits are only marginally weaker (despite S and T only being loop-induced by 
running). 

Our main focus was on the fermionic completion. We constructed an EFT valid in the 
limit V < rriy^ Mp and examined the 95% exclusion reach on the mass parameter space, 
for a variety of projected sensitivities for the precision electroweak and Higgsstrahlung 
measurements. Provided the coupling constant k is 0(1), we found that the precision 
electroweak limits are very powerful and primarily driven by T ~ {aK,/ae)S ^ S, owing 
to the violation of custodial symmetry. Already with current limits, the 95% confidence 
exclusion reach on Mp for <C Mp, is up to Mp ~ 675 GeV (2.9 TeV) for k = 1.0 (2.0). 
For the most optimistic projections we consider for the various possible configurations of 
the ILC, CEPC and FCC-ee colliders, these lower limits rise to 1.2 TeV (5.3 TeV), 1.8 TeV 
(7.7 TeV), and 2.0 TeV (8.5 TeV), respectively (see figure 2). 

For this model, the precision electroweak limits are generically more powerful than the 
precision Higgsstrahlung cross-section measurement due to the parametric enhancement of 
T over {AozhM^zh by 1/ae- Nevertheless, precise measurement of the latter also yields 
good complementary exclusion reach, which is useful because closely-allied UV completions 
in which the T parameter can be made small will therefore still be fairly strongly constrained 
by this measurement (other probes, such as /i —>■ 77 , may also become important). The most 
optimistic 95% confidence exclusion lower limits from Higgsstrahlung on Mj?, in the limit 
rriy^ <C Mp^ for the various possible configurations of the ILC, CEPC and FCC-ee colliders, 
are 490 GeV (1.5 TeV), 590 GeV (1.9 TeV), and 620 GeV (2.1 TeV), for k = 1.0 (2.0), 
respectively (see figure 3). 

We also considered the full one-loop computation of the electroweak precision observ¬ 
ables in the more general mass parameter space where this model does not necessarily 
provide a UV completion for the CP-even fermionic Higgs portal as we have defined it. The 
one-loop computation was found to agree well with the EFT computation for Mp. 

The model is significantly more constrained for the ‘opposite’ mass hierarchy rriy- ^ Mp 
owing to the fact that S and T are comparably large in this region; for example, we find 
95% confidence lower limits on for Mp = 100 GeV are typically a factor of 3-5 times 
higher than the corresponding lower limits on Mp for = 100 GeV (see figure 4). 

Additionally for the fermionic completion, the appearance in the EFT analysis of oper¬ 
ators modifying the coupling of the Higgs to SM fermions (see table 5) raises the prospect 
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of a modification to the rate for h ^ bb which will be measured to fairly high accuracy at 
lepton colliders; we defer investigation of this point to future work. 

Overall, we see that the sensitivities possible at future e'’'e“ machines through mea¬ 
surements of precision electroweak observables and the Higgsstrahlung cross-section will 
allow significant improvements in the exclusion reach for the CP-even fermionic Higgs por¬ 
tal over current limits, pushing into the (multi-)TeV range of particles masses, well beyond 
the direct reach of the LHC. 
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A One-loop computation of S and T 

The one-loop computation of S and T was performed in dimensional regularization in the 
MS scheme. The IJ self-energy (where I and J are SM vector bosons) arises from diagrams 
such as those in figure 5. For particles A and B with masses Ma and Mb, respectively, 
running in the loop, we find a contribution 

Ao(M2) + Ao(M|) - (M1 + M2) + _^2 Bo(p2,M2,M|) 

+^Bo(p2, Ml, Ml) {{Ma - Mb? - AMaMb) 

^ 1 Bo(pMu1,M|)-Bo(0,M1,M|) 

”^2 \ A B) 

-|- counter-terms 

(A.l) 

where Aq and Bq are the usual Passarino-Veltman scalar integrals [74]. Expanding n(p^) 
in powers of the gauge-boson momentum as n(p^) = n(0) + n'(0) , we find, at the 


n7j;AB(p^) = 


9iab ■ gjAB 


Figure 5. A one-loop contribution to the IJ self-energy function arising from particles A and B 
running in the loop. 


k + p 
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MS scale //, 


T-r giAB ■ gjAB 

n„;4B(0) = 


( M\- AMaMb + Ml 


+ 3 (M\ - Mi) 


0 3jS —y MQ , 


-1 


2MlMl log {Ml/Ml) 

-2Ml {6 MaMb - 3M2 + Ml) log (^7^1) 
+2Ml {QMaMb - 3Ml + Ml) log (/^VmI) 


/ 

(A.2) 

(A.3) 


and 

n/j;AB(0) = 


ff/AB • gjAB 


/ 

(Ml - M|) 




/ 6Ml(Ml-M|)log77M2) 

+ 6 M| (M| - M|) log [^‘^/Ml) 

\ +9M|Ms - 16M|M2 + 9MaMI + 2Mj + 2Ml 
\ +6MlMl {-SMaMb + Ml + Ml) log {mUmI) 


giAB • gjAB 
Idvr^ 


X ( -- ) log (/u7^7 as Ma^ Mb = M. 


/ 

(A.4) 

(A.5) 


In the broken phase of the theory, the mass-eigenstate basis of new fermions is {ni, n 2 , C}, 
where 


/ ni ' 

V^7 

Here the mixing angle is given by 


cos 6 — sin 0\ ( X 
sin 9 cos 0 / i A^, 


tan 20 = 


\/ 2 , 


KV 


Mp — 


and these particles have eigenmasses 


1 



1/2- 

Mi = - 

-|- Mp — 

{Mp — rriylf' + 2 k^v‘^ 


1 



1/2- 

M2 = - 
2 

1 

>< 

+ 

+ 

{Mp — m^)'^ + 2 k^v‘^ 



(A. 6 ) 


(A.7) 


, and Me = Mp. 


(A. 8 ) 

(A.9) 


The W~^W , ZZ, Zj and 77 self-energy functions and their derivatives can be immediately 
found by summing the contributions from eqs. (A.2)-(A.5) over all the allowed new-physics 
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Table 7. The new-physics particles {A,B) which run in the self-energy loop shown in figure 5. 


Self-energy IJ 

Particles (A, B) in loop 

w+w- 

(C,2) 

zz 

(1,1), (2,2), (1,2), (2,1), (C,C) 

Zl 

iC,C) 

77 

(C,C) 


particles which run in the loop (see table 7), using the following results for the couplings: 


gzcc = 


9Z11 = 


9^ - jg'? 

2v'FTFF’ 

Va^+igT . 2n 
- -^-sm e, 


9'yCC = 


gg 




Va‘^ + iaT 2 q 

9Z22 = -- --cos e, 


- a ■ n 

aw+ci — aw-ic — sm6, 

-an 

aw+c2 — aw-2c — 


J y^a^ + {a'f . a a 

and gzi2 = <7^21 = -smt/costl. 


(A. 10) 
(A.ll) 
(A.12) 
(A. 13) 


S and T can then be computed from their dehnitions, eqs. (4.4) and (4.9), and the 
results shown here; we have checked explicitly that all terms depending on g cancel in the 
results for S and T, as they should. Note that T arises as a small difference between two 
numerically large terms, and so care is required to obtain a numerically accurate result if 
these terms are not manually cancelled. 


B A clarification on matching in the presence of mixing 

This Appendix aims to clarify a point about the matching at the scale g = Mp which 
may not be immediately familiar to the reader: when {H) ^ 0, the mass-eigenstate y held 
we write in the intermediate-energy EFT^^ should not be thought to be identihed with 
the (mass-eigenstate admixture) y held in the UV theory; rather the x held in the EFT 
should instead be identihed with the light-mass-eigenstate ni in the UV theory. This fact 
notwithstanding, it is common practice, although something of an abuse of notation, to still 
call the held in the EFT ‘x’- 

Consider hrst a general argument. We have presented both the EFT and UV theory 
in terms of SU{2) x ?7(l)y-symmetric operators. However, because (a) the limit {H) —)• 0 
is smooth in the sense that the F doublet in the UV theory remains sufficiently massive 
{Mp ^ Kv, m^) to be integrated out, and that the EFT power-counting is in p{H)/Mp, so 
no EFT operators or Wilson coefficients diverge as {H) —?■ 0; (b) the only source of EWSB 
in either the UV theory or the EFT is the existence of a non-zero (77); and (c) the EFT 
is explicitly constructed to match the amplitudes of the UV theory at fi = Mp, it follows 

this Appendix, ‘EFT’ always refers to the intermediate-energy EFT, after the F field is integrated 
out at /i = Mf, but before the x held is integrated out a.t jj, = m^. This is the EFT described by the 
operators and coefficients in Tables 2 and 3. 
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that when (H) ^ 0, the EFT must by construction capture all the (ff)-dependent effects 
which are present in the UV theory, at each order in the EFT power counting. 

Put another way, up to the truncation order of the EFT, one obtains the same physical 
predictions from the EFT by (a) matching in the broken phase {{H) 7 ^ 0) of the UV-theory 
by integrating out the heavy-mass-eigenstates n 2 and C, and writing the EFT in terms 
of the light-eigenstate ni and the physical Higgs field h] or (b) matching in the unbroken 
phase {{H) = 0), integrating out the F doublet whole and writing the EFT in terms of 
(the EFT field) y and the S'C/(2)-doublet field H, and then turning on a non-zero {H) in 
the EFT after the matching. Once {H) is turned on, the EFT automatically includes the 
corrections to the x field which allow its identification with the same physical degree of 
freedom represented by the light-mass-eigenstate ni of the UV theory. 

We emphasize that it is by construction that the EFT must reproduce the UV-theory 
amplitudes and so must capture the mixing of the fermions x N in the UV-theory. 
In the S'C/(2)-symmetric EFT, this mixing is present in the form of the higher-dimension 
operators of the EFT. We now show examples of this by explicit computation. 


B.l Example 1: two-point functions of the UV-theory ni field and the EFT 
field X 

The clearest demonstration that the EFT field x represents the same physical degree of 
freedom as the UV-theory field ni is to show that these two fields have the same two- 
point function; or, equivalently, that when they both have canonical kinetic terms in their 
respective Lagrangians, they have the same physical mass (up to corrections at the truncated 
order in the EFT power-counting). 

The UV-theory ni field has a canonical kinetic term and mass given in eq. (A. 8 ). 
Expanding this result out for Mj? S> Kv,my^ gives 


Ml 




2Mp 


V 

2 M2 




+ 


4M|, 


+ 0{Mp 


-4i 


(B.l) 


The EFT x field does not have canonical kinetic terms. The relevant higher-dimensional 
operators in the EFT are (see Tables 2 and 3) 


/:eft V ix$x - m^xx + ^ Oes + ^ O7A + ©(M"") (B.2) 

r J? p' 

^2 1 ^2 


CQB 


C7A 


= ix0X - m,^XX + ~ ^-C-) 


' 2 Ml 


1 M 


- + h.c.) + 0 (M/). 

F 


(B.3) 


Breaking electroweak symmetry by setting {H) = (0,r;/\/2)^, the relevant terms become 


£eft V i-X^X 






^ — —V 

' 2M2_ 

- XX 

2Mp 

AMl 

XDX+XDX 


+ 0{M 




(B.4) 
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This can be re-cast to canonical form by performing a field-redefinition, which does not 
change the physical content of the theory: 


X 


1 - 


4M| 




4M|, J 


i$X 


K V l ^ 


(B.5) 


In terms of the redefined held, we have 
-Ceft ^ 'ix$X - XX 


m-yK^'v 


2Mp 




2Mj, 


4M|. 


+ 0(M/). (B.6) 


Comparing eqs. (B.l) and (B.6), we have established the desired result. 


B.2 Example 2: Z-boson coupling to the EFT field x 

We now look at an example for the EFT x interaction terms. Consider that the UV-theory 
X held is an exact SM-singlet, while N G F has SM gauge-couplings. When (FI) ^ 0, 
the mass-basis rotation/mixing, eq. (A.6), generates SM gauge-couplings for both the ni ,2 
mass-eigenstate fermions; see eqs. (A.10) to (A.13). 

If the EFT held x were to be identihed exactly with the x held in the UV theory, it 
would necessarily follow that the former should have no SM gauge-couplings; if instead the 
EFT held x were identihed with the light-mass-eigenstate ni of the UV theory, we would 
expect the former to couple to, e.g., the Z-boson. The second alternative is the correct one, 
as the following short argument shows. 

Consider that in the UV theory the coupling constant of the ni vector current to the 
Z-boson is gzii = 6 where tan20 = x/2kv/{M p — see eqs. (A.7) 

and (A. 11). In the limit where Mp ^ kv, nix, we have 

9ZU ~ + {a'f + • • • ■ (B.7) 

Meanwhile, in the EFT the following operator is present: 

/:eft D ^ O,a D -\V9^ + W ^ X7/.X + • • • • (B.8) 

(Note: the field-redehnition in Example 1 above does not modify this leading-order cou¬ 
pling.) One can readily see that the leading-order coupling of the Z-boson with the vector- 
current of the EFT x held has the same sign and magnitude, and appears at the same 
order in the EFT power-counting, as the leading-order coupling of the Z-boson with the 
vector-current of the UV-theory m held. Again, we draw the reader’s attention to the fact 
that it is a higher-dimension operator in the EFT that has captured an effect of UV-theory 
fermion mixing when {H) ^ 0. 
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